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We provide a new approach, along with extensions, to results 
in two important papers of Worsley, Siegmund and coworkers closely 
tied to the statistical analysis of fMRI (functional magnetic resonance 
imaging) brain data. These papers studied approximations for the ex- 
ceedence probabilities of scale and rotation space random fields, the 
latter playing an important role in the statistical analysis of fMRI 
data. The techniques used there came either from the Euler char- 
acteristic heuristic or via tube formulae, and to a large extent were 
carefully attuned to the specific examples of the paper. 

This paper treats the same problem, but via calculations aris- 
ing from the so-called Gaussian kinematic formula. This allows for 
extensions of the Worsley-Siegmund results to a wide class of non- 
Gaussian cases. In addition, it allows one to obtain results for rotation 
space random fields in any dimension via reasonably straightforward 
Riemannian geometric calculations, whereas previously only the two- 
dimensional case could be covered, and then only via computer alge- 
bra. 

Finally, by adopting this more structured approach to this par- 
ticular problem, a solution path for other, related problems becomes 
clearer. 

1. Introduction. In this paper we provide a new approach, along with 
extensions, to results in two important papers closely tied to the statistical 
analysis of fMRI (functional magnetic resonance imaging) brain data. The 
first, by David Siegmund and the late Keith Worsley [9], treated the problem 
of testing for a spatial signal with unknown location and scale perturbed by 
a Gaussian random field on M N . While their motivation came from fMRI, 
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their results are also relevant to areas such as astronomy and genetics, or 
indeed to many applications in which one searches a large space of noisy 
data for a relatively small number of signals which manifest as 'bumps' in 
the random noise. 

In a second paper, this time together with Shane and Sigal [8], they treated 
an extended version of the same basic problem, this time allowing the asym- 
metric bumps to also have a direction associated to them, in addition to 
scale changes that are direction dependent. This added an additional 'rota- 
tion parameter' to the model, and complicated it to the extent that they were 
unable to treat anything beyond the two-dimensional (JV = 2) case. Thus, 
for example, while their results could be used to examine cortical surface 
images, they could not be applied to three dimensional fMRI analysis. 

For obvious etymological reasons, which will become mathematically clear- 
er soon, the first of these cases was called the scale space case, and the second 
the rotation space case. In both cases, the mathematics of these papers con- 
centrated on approximating the exceedence probabilities 



for large u, a smooth, Gaussian, real valued, random field /, and some 
parameter space M. For reasons that will become clearer soon, in the N- 
dimensional scale space case, the effective dimension of M was N + 1, while 
in the rotation space case it jumped to N + N(N + l)/2 (or 5, for the only 
case they could handle, when N = 2). 

The aims of the current paper are two- fold. Firstly, we shall generalize 
the results of the rotation space to any dimension, but by a very different 
method. Rather than taking a 'from first principles' approach, as was the 
case in [8], we shall compute everything from a basic result known as the 
Gaussian kinematic formula (GKF, cf. (2.18) below), using techniques of 
Riemannian geometry. In fact, we shall derive the rotation space results from 
a more general setting, which is of interest in itself. This is quite different 
from the analysis in [8], in which calculations were so involved that they 
could only be handled by symbolic calculus using computers, and, then 
as already mentioned, only in two dimensions. According to the 'no free 
lunch' principle, there is a price to pay here, in that the specific N = 2 
approximations of [8] have a slightly higher level of precision. 

Secondly, we shall recover the results of the scale space case, again using 
the GKF. We shall do this in two ways. Firstly, we shall obtain it as an 
immediate corollary to the 'more general setting' mentioned above. This is 
quick, but once again does not give the best results. However, we shall also 
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derive the 'best' results from a more detailed calculation, something which 
seems too difficult to do for rotation space fields. 

At first this might seem as if we are attempting to kill the proverbial ant 
with a sledge hammer, but it turns out that there are at least three advan- 
tages to this approach. First of all, the ant in question is a rather large one, 
and any way of dealing with it requires considerable effort. The basic calcu- 
lations in both [9] and here take about 10 pages, and the essential difference 
is whether one prefers to work with Riemannian topology or integral geom- 
etry, hard calculus, and occassionly symbolic (machine) calculus. This may 
be more a question of mathematical taste than anything else. The second 
advantage of using the GKF as the basic starting point, which is significant 
and no longer merely an issue of taste, is that while the results of [8] and [9] 
apply, as stated, only to Gaussian fields, the GKF allows immediate exten- 
sion to a much wider class of examples. We shall describe this briefly below 
in Section 2.5. 

Finally, both of these cases are instructive in terms of how to use the GKF 
in more general settings, and our hope is that presenting them in detail will 
provide examples which will make the GKF easier to use in the future for 
other, perhaps quite different, examples. 

We shall say no more in the way of motivation for this paper. Discussions 
and examples can be found in both [8] and [9] as well as the papers cited 
there, and the many more recent papers citing [8] and [9]. In particular, it is 
explained there why (1.1) is an important quantity in statistical hypothesis 
testing. Rather, we shall, in the following section, immediately turn to formal 
definitions of scale and rotation random fields, a description of the GKF and 
its relation to [8] and [9] , and a statement of the main results of the paper 
in Section 2.6. Fuller results for scale space fields, which go considerably 
beyond those of Section 2.6, are given at the end of the paper, in Section 8. 

Understanding and using any of the results of the paper do not require 
a deeper knowledge of differential geometry - the main tool of the proofs 
- beyond that which is developed prior to them being given. This is not 
true of the proofs, and in Section 4 we point to what else is needed. The 
remainder of the paper is taken up with proofs. 

2. Random fields, the GKF, and our main results. The common 
element behind all the random fields in this paper is Gaussian white noise. 

2.1. Gaussian white noise. The simplest way to define Gaussian white 
noise W on is as a random field indexed by the Borel sets B , satis- 
fying the properties that, for all A, B G B N with finite Lebesgue measures 
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X N (A) = \A\ and X N (B) = \B\, 

(2.1) W(A) ~ N(0, \A\). 

(2.2) An B = ® ^ W(A(J B) = W(A) + W(B) a.s. 

(2.3) vlnB = N and W(-B) are independent. 

It is then standard fare to define integrals 

(2.4) / <p(t)W(dt), 

Jm. n 



for deterministic ip £ L2(M ), and to check, as a consequence of (2.1)— (2.3), 
that such integrals are mean zero, Gaussian, and with covariances given by 

(2.5) e{ f <p(t)W(dt) [ ip(t)W(dt)\ = [ <p(t)il>(t)v(dt). 

[Jr n Jr n J Jr n 

2.2. Scale space fields. The Gaussian scale space random field is obtained 
by smoothing white noise with a spatial filter over a range of filter widths 
or scales. Formally, for a 'filter' h £ Li2(R N ) it is defined as 

(2.6) f(a, t) = a' N / 2 [ h (t^i) dW (u). 

JR." V a J 

Typically t ranges over a nice subset T of M. N , and we take a G [ai, a u ] , for 
finite < a i < o u . 

The parameter space therefore includes both scale and location parame- 
ters [9, 15]. From a statistical point of view, the scale space field is a con- 
tinuous wavelet transform of white noise that is designed to be powerful at 
detecting a localised signal of unknown spatial scale and location. 

This concludes the definition, and we now turn to some additional as- 
sumptions, as well as setting up some notation and further definitions that 
we shall require to state the main results and to carry out the proofs. 

For convenience, and without loss of generality, we shall assume that the 
filter h is normalised so that J RN h 2 {u)du = 1, which will ensure that the 
variance of / is also one. We shall also assume that h is twice continuously 
differentiable on R , and that all partial derivatives are in L2OR )■ 

With some loss of generality, that, in practice, is not generally too restric- 
tive, we shall assume that there exists a A > such that 



(2.7) 



/ Vh(u)(Vh(u))' du = XInxN- 
Jm. n 
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(Note: We take all our vectors to be column vectors.) Equation (2.7) will 
follow, for example, if the filter is spherically symmetric. Two common ex- 
amples are given by the standard Gaussian kernel 



1 



"W = ^- |,|3/2 . 



and the Marr or 'Mexican hat' wavelet 

AN 



h{t) 



(N + 2)ir N / 2 



1/2 / \ f \2 



-1*172 



It is easy to check that in the Gaussian case (2.7) holds with A = 1/2 and in 
the Marr wavelet case with A = (N + 4)/(2N). Furthermore, the covariance 
function of the scale space field is easily calculated, via (2.5), to be 

C^MU^M)) = r^rwrJ h ( f -^) h (**^) du 

{cricr 2 r /2 Jrn \ ai J \ <J 2 J 

1 f h( tl - h - V )h(^\dv. 



{aia 2 ) N / 2 Jrn V °\ ) V °2 

Note that, for fixed a, it follows from the last line above that / is stationary 
in t. However, / is most definitely not stationary as a process in the pair 
(a,t). 

To apply the GKF, and to state our results, we shall need information 
about the derivatives of /, which exist due to our assumptions on h. Either 
from first principles, by differentiating the covariance function, or by using 
(5.5.5) of [3], we find that all first order partial derivatives of / with respect 
to the space variables are uncorrelated with the first order derivative in 
the scale variable. This lack of correlation, which of course is equivalent 
to independence since all fields are Gaussian, is crucial to the proofs, and 
without it it is unlikely that we could complete the detailed calculations that 
we shall carry out later. The same techniques also give 

(2.8) k = a 2 ■ Var ( = f [(u,Vh(u)) + Nh{u)/2] 2 du, 

\ da J JnN 



(2.9) K = Var(V t /(<7,t)) = a 2 [ Vh(u)(Vh(u))' 



du, 



where Vt/ denotes the gradient of / with respect to the elements of t only. In 
view of (2.7), we have A CT = <j 2 \In x n- We should note that it is not obvious 
at this point that the quantity a 2 Var (df (a, t)/da) should be independent 
of a. However, this is in fact a consequence of calculations in Section 6. For 
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the Gaussian and Marr wavelet kernels the values of k are 1/2 and N/2, 
respectively. 

Excluding Section 8, from now on we shall assume that h is spherically 
symmetric, so that we can define k : M+ — > K by 

(2.10) h(x) =k(\\x\\ 2 ). 

In this case, A, k and other constants which will appear throughout can be 
expressed in terms of k. In particular, retaining consistency with (2.7) and 
(2.8) and introducing three new terms we have 



^2,2 



(2.11) 



A = 4 / k 2 (u'u) u\ du 

I ' 2 / / \ 2 2 

/ k (uu) UiU 2 du 
Jr n 

«4 = / k (uu) uf du 



and also define 

(2-12) C\ N k K - = N 2 (1 + P- + K2)2 ) + 2N K2 , 2 . 

Note for later use that 

(2.13) K4 — 3^2,2 = / k 2 (u'u) (uf — Zulu 2 ,) du = 0. 

Jr n 

This can be easily proved for a general isotropic kernel (instead of k 2 (u'u) 
above) by using Fubini's theorem to restrict to the 2-dimensional case and 
then expressing the integral using spherical coordinates. 

2.3. Rotation space fields. As already described, rotation space random 
fields are based on the same underlying white noise model as scale space 
fields, but also allow for rotations and scale changes that are direction de- 
pendent. More precisely, write 

S(a 2 2 , a~ 2 ) 4 {S g Sym+ : a 2 2 I <S< a~ 2 l} 

where Symjy C Sym^ is the set of N x N positive definite symmetric matri- 
ces, < <7i < 0- 2 < oo, K = N(N + l)/2, D = N + K and for two matrices 
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A and B we write A < B if, and only if, B — A is non-negative definite. 
Then the rotation space random field, defined over the D-dimensional space 
M = T x §, is given by 

(2.14) f(t,S) = \S\* [ h(s^ (t-u)) dW (u) . 

Jr n v 1 

Again, we shall assume that h(t) is spherically symmetric, twice continu- 
ously differentiable on K , that all partial derivatives are in L 2 (R N ). and 
normalised so that J^ N h 2 (u) du = 1. Note that the eigenvalues and eigenvec- 
tors of S determine the scalings and the corresponding 'principal directions' 
for the rotation space field. 

The constant analogous to Ck,N m rotation space is 

(2.15) C k , N ± (2^- N Up^c k , N , 
where 

(2.16) c\ N t (2k 2 , 2 ) K (l + N {2^ 2 y 1 + f + «2, 2 

Reassuringly, when iV = 1, in which case scale space is equivalent to 
rotation space, it is easy to check that Ck = Ck- Another constant which 
will appear later is 

(2-17) D% N t 1 + | + 3k 2 , 2 . 

As for scale space fields, we shall once again need information about 
variances and covariances of the first order derivatives of rotation space 
fields. This time, however, the formulae are somewhat more complicated, 
and, since one of the parameters of f(t,S) is a matrix, we shall have to 
explain what we mean by differentiation with respect to this parameter. 
Details will be given below in Section 5.1. At this point we note only that, 
as before, derivatives with respect to t and S are independent, and, again, 
this is crucial to our ability to carry out detailed computations. Furthermore, 
although once again / is neither stationary nor isotropic, it does have zero 
mean and constant variance. 

2.4. Regularity conditions. Before we state the GKF, we need to say 
something about the regularity conditions we need for it to hold. These are of 
three kinds, relating to the smoothness of the scale and rotation space fields, 
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the structure of T (and so of M), and the smoothness of a transformation F 
that we shall meet only in the next subsection and so defer until then. Full 
details are given in Chapters 11 and 12 of [3], and you will have to turn to 
there to see why the conditions of this subsection suffice for the GKF of the 
following one. 

We start with T, which we take to be a compact domain 1 in M. with a C 2 
boundary. For the remainder of this paper we require no convexity conditions 
on T, but we do point out that in order to apply the Euler characteristic 
heuristic of Section 3 it is necessary to assume that T is locally convex in 
the sense of Definition 8.2.1 of [3]. 

The second set of conditions are required to ensure that / is, with prob- 
ability one, C 2 , as a function of each of the space and scale or rotation 
variables, and also non-degenerate in a certain (Morse theoretic) sense. Suf- 
ficient conditions under which this holds are given in [3]. It is easy to check 
that these conditions are satisfied for the Gaussian, Marr, and other smooth 
kernels. 

2.5. The Gaussian kinematic formula. We shall state the GKF in more 
generality than we need it for this paper, since it requires little more than 
an additional sentence or two to do so. Suppose that M is a C 2 , Whitney 
stratified manifold satisfying some mild side conditions and D a similarly 
nice stratified submanifold of M k . For the exact definitions of 'mild' and 'nice' 
see [3], but, for the purposes of this paper, it suffices that M is either a scale 
or rotation space, with T satisfying the conditions of the previous subsection. 
Also, if F : R k -»• R is C 2 , then it will suffice for us that D = F~ l ([u, oo)) 
for some u. 

Let / = (f 1 , . . . , f k ) : M — > M fc be a vector valued random process, the 
components of which are independent, identically distributed, real valued, 
C 2 , centered, unit variance, Gaussian processes, satisfying the conditions of 
the previous subsection. The Gaussian kinematic formula (due originally to 
Taylor in [10, 11] and extended in [3, 12]) states that, for certain additive 
set functionals Co, ... , £dimM, 

dim M—i 

(2.18)E{£ l (Mnr 1 (D))}= E 

j=0 

x We could actually assume less, but we avoid doing so to keep the proofs of the paper 
relatively transparent. In particular, it would suffice to assume that T is a regular stratified 
manifold in the sense of Definition 9.2.2 of [3], viz. T is a C 2 Whitney stratified manifold 
satisfying some additional, mild, regularity conditions. 



{2^l 2 C l+J {M)M k {D). 
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There is a lot to explain here. The easiest are the combinatorial 'flag' coef- 
ficients 



j, 



where oj n is the volume of the unit ball in W 1 . The A4j(D), known as the 
Gaussian Minkowski functionals of D, are determined via the tube expansion 

(2.19) p{iV(0,4 xfc )G{yG]R fe : \\y, D\\ < p}} = £ ^M){D). 

In the case D = F~ 1 ([ti, oo)) it typically involves no more than multivariate 
calculus to compute the Mj(D), and many examples are given in [3] and 
[4]. We shall give one example in a moment. 

The most important element of (2.18) are the Cj, known as Lipschitz- 
Killing curvatures (LKCs). The first of these is rather simple, as it is always 
the Euler characteristic x; i-e. 

C (A) = X (A) 

for all nice A. Thus, (2.18), with i = 0, provides us with an expression for 
E{x (Mn/ _1 (D))}, which is what is needed to approximate exceedence 
probabilities, as explained in the following section. 

For a simple but extremely important example, suppose that k = 1 and 
F is the identity function. Then it is easy to compute the .MjQii, oo)) to 
see that, with i = 0, and / : M — > M. satisfying the conditions above, (2.18) 
becomes 

dimM 

(2.20) E{ X (A u (f,M))} = e-^ (2^)"^ Cj (M) (x) . 

i=o 

Here 

(2.21) A u (f,M) ={teM: /(*)>«} 

is an excursion set, the H n , for n > 0, are Hermite polynomials, and 
H-\(x) = [(^(x)] -1 tp(u) du, where 93 is the standard Gaussian density. 

Thus, we shall have everything we need, once we know how to compute 
the Lipschitz-Killing curvatures, which is what this paper is about, for scale 
and rotation space fields 2 . 



2 Note the important fact that the structure of M affects only the LKCs, but neither 
the Minkowski functionals nor the flag coefficients. 
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First of all, a word to the purely applied statistician who would prefer to 
leave us now and turn to a numerical, data driven method for estimating 
LKCs from data. Two different methods [2, 14] for doing this are currently 
available. For the remaining readers, we now turn to a more technical de- 
scription of LKCs. 

For this, we assume that M is a regular (in the sense of [3]) stratified 
manifold. These are basically sets that can be partitioned into a disjoint 
union of C 2 manifolds, so that we can write 

dim M 

(2.22) M = | djM, 

3=0 

where each nonempty stratum, djM, < j < dim(M), is a j-dimensional 
manifold whose closure contains d{M for all i < j. For a iV-dimensional 
cube, d]yM is its interior, d n _\M is the union of its (N — l)-dimensional 
(open) faces, and so on, down to d$M, which is the collection of its vertices. 
We shall see how this works for scale and rotation space in the following 
sections. 

If such a space is locally convex, then a simple way to define (a Euclidean 
version of) the LKCs is via a tube formula, which states that, for small 
enough p, and M C E , 

(2.23) U E N (\^t G R N : inJ \\t - s\\ < p^ = Yj^N-jp N ~ j ^f (M), 

where (for reasons that will become clearer later) we write %fj for Lebesgue 
(Hausdorff) measure on Mr . From this it is easy to check, for example, that, 
for rectangles, ^(H^ilO, Ti]) = ■ ■ -Ty, where the sum is taken over 

the ( ■ ) different choices of subscripts i±, . . . Similarly, for A^-dimensional 
balls of radius A, Cf(B%) = (^)A J W/^_j-. 

Note that it will always be true that Cfj(M) is the volume of M, while 
Cfj_ l (M) is half the surface measure of M. Furthermore, Cq(M) = x(M). 

In what follows, however, we shall need to work in a non-Euclidean set- 
ting, as the Lipschitz-Killing curvatures in the GKF, and thus in (2.20), are 
computed relative to a Riemannian metric induced by a Gaussian field. This 
metric g on M is defined by 

(2.24) g x {X x ,Y x ) i E{(X x f x )(Y x f x )}, 

for any i and for X X ,Y X G T X M, the tangent space to M at x G M. This 
can also be written in terms of the covariance function as 

(2.25) g x (X X ,Y X ) = X x Y y C(x,y)\. 
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Both of the Gaussian fields of interest to us - viz. the scale and rotation 
space fields - which we defined on compact sets, can easily be extended, to 
1^ x R and x Sym^, respectively, and so (2.24) defines metrics on these 
spaces. In all that follows, unless stated otherwise, these metrics, or their 
restrictions, are what we shall work with, and it will be in relation to them 
that Hausdorff measures T-Lj and Lipschitz-Killing curvatures Cj are defined. 
When we consider a submanifold L of either of these spaces the notation 
Hj, £j will be used. When the standard Euclidean metric is used these will 
be replaced by T~L^ and £,? . 

A Riemannian metric defines a volume form, a Riemannian curvature, 
and a second fundamental form on M (and, indeed, on each the djM). 
Then CdimM(M) is its Riemannian volume, and so is an integral of the 
volume form. £dimM-i(-^) is closely related to the (Riemannian) area of 
^dimM-i(-^O) the (dimM — l)-dimensional boundary of M. The remaining 
Lipschitz-Killing curvatures involve integrals of both curvatures and second 
fundamental forms, and are more complicated. Rather than give the general 
form (for which you can see Chapter 10 of [3]) we shall give specific defini- 
tions for each of the cases we treat, in the following sections. Now, however, 
we shall give (relatively) simple expressions for the cases of interest to us. 
These are the main results of this paper. 

2.6. Main results. We shall present the results of this section in stages. 
In the first we shall represent the Lipschitz-Killing curvatures as an integral 
for a situation more general than that needed to handle rotation space fields. 
This is Theorem 2.1, and it shows how to write these integrals as products 
of Lebesgue measures of subsets of M. N and integrals over submanifolds of 
Symi. In Theorem 2.2 we shall compute these integrals for certain situa- 
tions, which will lead to Corollaries 2.3 and 2.4. These give relatively simple, 
quite explicit, expressions for the highest order Lipschitz-Killing curvatures 
for rotation and scale space fields, respectively. More detailed information 
for scale space fields is given in Section 8. 

We say that a manifold L C Sym^ is orthogonally invariant if S £ L => 
QSQ' G L for every Q £ O(N), the set of orthogonal matrices. 

For such L of dimension to, we define 



(2.26) 




where 



0{N) 



(QSQ') jxj 2 du (Q) . 
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Here u> is the Haar probability measure on O (N) and for a matrix A, Aj X j 
denotes the corresponding upper left submatrix of size j x j. The measure 
is the Haar measure induced by the restriction of the metric induced by 
the field g to a fiber t x Syihjy, which, as we shall see, is independent of the 
choice of t. 

With these definitions, our first result is the following. 

Theorem 2.1. Let f : A[ x B k C 1^ x Syrri^ — > R be the restriction of 
the rotation space random field defined in Section 2.3, to the product of an l- 
dimensional regular stratified manifold A[ C R N and a k- dimensional regular 
stratified manifold C Syrn^. Furthermore, suppose that each stratum 
djBk, < j < k, is orthogonally invariant, and that the regularity conditions 
of Section 2.4 vis a vis f all hold. Then, 

C l+k {A t x B k ) = X^nfidiAt) F kil (d k B k ), 

and 

A +fe _i(M) = l^-^nt^d^A^F^dkBu) 

+ X^nfid^F^idk^Bk)], 

where A is given by (2.7) or (2.11). 

It is not hard to see that TZj (S) is a symmetric function of the eigenvalues 
of 5. Furthermore, as we shall see, as each stratum of B m is orthogonally in- 
variant, its Hausdorff measure is in fact invariant under S \— > QSQ 1 , for any 
orthogonal Q. Hence, the F m j are expressible as integrals over certain func- 
tions of the eigenvalues of S £ d m B m . In particular, for rotation space, with 
search region T x 8(0>7 2 , cf 2 ), we arrive at the following explicit expressions. 

Theorem 2.2. With the notation defined above, under the conditions of 
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Theorem 2.1, and for a spherically symmetric kernel, 

r Jl 

F KtN (§(af, err 2 )) =C k , N 2 | A (A)| J] A^ 2 dX 

J a 2 <X\<---<Xm<o- 1 j—i 

F^.^S^" 2 ,^ 2 )) = C KN I _ 2 |A(A)|ftjv-i (diag(X)) 

^ „ Af-1 

^_ 1>a ksk 2 ,o) = §^ / 2 2 ia(a)i n ^ m 



'fc,Af J<t 2 "^<A 1 <---<A jv _ 1 <o- ] 7 

Af-1 

_fV-2 TT U -21 , JV-2 TT 1-2 \ I 

2 lll A i- ff 2 1+^1 11 Tl - X 3\ 

3=1 3=1 



dX 



where Ck t N an d D^^n are given in (2.15) and (2.17). 

Theorem 2.1 implies the following corollary, with all F m j's as given in 
Theorem 2.2. 

Corollary 2.3. Under the conditions of Theorem 2.2 and the regularity 
conditions of Section 2.4, the top two Lipschitz-Killing curvatures for the 
rotation random field on T x §(o"^ 2 ,o"f 2 ) are 

£ D (Tx% 2 V[ 2 )) = C N+K (T x§(af,a^ 2 )) 

= X N ^ 2 n§(d N T)F KN (d K S(a^,a^)) 

and 

£;,_i(Tx% 2 Vi 2 )) 



l^ A (Af-l)/2^_ i(a7v iT)jP ^_ i( ^g ((J -2 )CJ -2 )) 

+ \ N ' 2 H% (d N T)F K . hN (d K ^S(a^ 2 , exf 2 )) 



To obtain a parallel result for scale space fields, we use the fact that the 
scale space random field on T x [a\, 0-2], with a spherically symmetric kernel 
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has the same distribution as the restriction of the rotation space random 
field defined in Section 2.3 to T x B(<r^~ 2 , crj" 2 ), where 

B(a2 2 ,a^ 2 ) = {ul : ovT 2 <v < a^ 2 } C S(a^ 2 , a^ 2 ). 

Once we study the metric induced by the rotation field, Theorem 2.1 will 
effortlessly yield the following result. 

Corollary 2.4. Under the conditions of Theorem 2.2 and the regularity 
conditions of Section 2.4, the top two Lipschitz-Killing curvatures for the 
scale space random field on T x [o"i,<72]> are 

C N+l {T x D^-Vr 2 )) = \ N/2 U%(d N T) F liJV (diDKVr 2 )) 



(2.27) 



C k , N \ N ' 2 U E N (d N T) / u N '^ du 



-2 

'1 



2 
2 



and 

C N (T xD(^ 2 l( rr 2 )) 
1 
2 



(2.28) = \{ A^-^f.i^i^Fi^x^BK 2 ,^ 2 )) 

+\ N / 2 n§(d N T)F , N (d B(a2 2 ,aY 2 )) 



{ -(N~l) _ -(N-l)\ 

= i(^-^s. 1( w 1 {N :^ /2 > 

where C^.n is as i- n (2-12). 

In Section 8, we analyse scale space fields independently of the rotation 
space case for a manifold with a boundary and without assuming spherical 
symmetry. The formulae for all Lipschitz-Killing curvatures are given via 
(8.2), the terms of which are all computed in Section 8.7. 

2.7. Comparison with existing results. While the formulae of Theorem 
2.2 may look rather forbidding, one should note two facts. The first is that, 
for dimension N > 2, they are the first of their kind. The second is that 
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although for dimension N = 2 they are implicit in [8], one computation re- 
lied on computer algebra and gave a large, unstructured formula, while the 
second, based on tube formulae, does not precisely give the LKCs. (This is 
not unreasonable, since that computation was directed explicitly at approx- 
imating exceedence probabilities, and not computing LKCs.) To recoup the 
results of [8], we treat, as an example, the case N = 2 (so that D = 5). In 
that case, it is possible to carry out the integrations needed to compute £5 
and £4. A page or so of calculus (irritating but standard) shows that, for 
the Gaussian kernel, 

(2.29) £ 5 (M) = 2" 4 7tcj 1 - 2 \T\ [r 2 - 1 - (r 2 + l) lnr] , 

where \T\ is the simple Euclidean area of T and r = a\/a2- 

The LKC £4 is a little harder to calculate, since while the first term inside 
the brackets is again amenable to simple calculus, the second term leads us 
to special functions. However, a simple answer is attainable, and 

£4 (M) = 2- 7/2 trt^ 2 \T\ (r 2 - 1) lnr 

+ 2" 9 / 2 7T \dT\ ( (A 2 - Ai) (AiA 2 )~ 3 / 2 K (Ai, A 2 ) dX, 

where \dT\ is the length of dT, and, with E the elliptic integral E (y) = 
£ /2 s/1 -y 2 sin 2 BdS, we have ft (Ai,A 2 ) = vrAi^ ((A 2 - A|)/A 2 ) /2. 

Note also that, as stated, Theorem 2.2 is weaker than the result for the 
scale case in Section 8 in that that it gives only the highest two LKCs, rather 
than all of them. However, while more would obviously be desirable, it is the 
higher ones that are most important for the Euler characteristic heuristic 
discussed in the following section. We shall explain there why this is the 
case. 

3. The Euler characteristic heuristic: ECH. As intimated in the 
Introduction, investigating a signal plus noise problem by testing the null 
hypothesis of no signal often boils down to being able to compute the ex- 
ceedence probabilities (1.1). When dealing with smooth Gaussian processes, 
there is typically no way to directly compute these probabilities, and they 
must be approximated. Generally, one requires an approximation that is 
good in the tail of the extremal distribution, where by 'good' we mean that 
in the last 5% or so of the tail the estimates of probabilities or threshold 
levels have relative errors of only a few percent. 

There are a number of ways to carry out such approximations, and two of 
these are at the center of the papers [8] and [9] which were discussed in the 
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Introduction. One of these relies on what are known as tube formulae, and 
moves the problem from one defined on Gaussian processes on Euclidean sets 
M to one involving computation of volumes of tubes around embeddings of 
M in high dimensional spheres. 

The Euler characteristic heuristic (cf. [1, 3]) provides what is generally 
an easier route to these approximations, and has the advantage that it also 
provides information on the sample paths of random fields at all levels. The 
ECH states that for smooth random fields /, on nice, locally convex, compact 
sets M, 

(3.1) p{ sup /(*)>«} « E{ X (A u (f,M))}, 
for large u. 

In the case of smooth, Gaussian, /, with zero mean and constant unit 
variance (but with no stationarity assumptions) the approximation in (3.1) 
can be quantified (cf. [3, 13]), and it is known that 

(3.2) liminf-^ 2 log|Diff /iM (u)| > \ U + ) , 

where we write Diffj m(^) f° r the difference of the two terms in (3.1), and 
a c(f) > is an /-dependent constant that is, in many cases, computable. 

The importance of (3.2) is that it shows that the two terms in (3.1), both 
of which are known from the general theory of Gaussian processes to be of 
order ii dimA/_1 exp(— tt 2 /2), differ by a term of exponentially smaller order; 
viz. of order exp(— u 2 /2(l + <r 2 (/)). This fact is what justifies the use of the 
ECH for Gaussian fields (in which case it is no longer heuristic) and in large 
part motivates the computations of this paper. 

It also explains the claim of the previous section, that the leading Lipschitz- 
Killing curvatures are the important ones in applying the ECH, since these 
are the coefficients of the highest order terms in the approximation. In the 
case of real valued Gaussian fields this is clear from the ECH, the expansion 
(2.20), and properties of Hermite polynomials. 

4. Before we start the proofs. As mentioned in the Introduction, it 
should have been possible to read the earlier sections of this paper based 
on what was defined there. No real knowledge of differential geometry was 
required. That will not be true for the proofs. 

The derivations of the leading Lipschitz-Killing curvatures for the rotation 
space case, which is what we have concentrated on so far, and all LKCs 
for the scale space case, which will be treated in Section 8, begin in the 
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same fashion. Each starts by taking the parameter space M, writing it as a 
stratified manifold, and carefully identifying the strata djM of (2.22). The 
next step will be to take the generic formulae for LKCs, given in (necessary, 
but sometimes painful) detail in [3] , and rewrite them for the current setting; 
i.e. for the specific stratification, and for the induced Riemannian metric 
(2.24). We shall not explain in detail how this is done other than to say, at 
each stage, which generic formula we are using. We leave it to the reader to 
verify that the generic formulae simplify to the ones we give here. Then, the 
Riemannian metric induced by the field is computed. 

At this point, however, the proofs diverge. For the scale space we compute 
all LKCs, and so geometric objects such as curvature and second fundamen- 
tal forms must be computed. For the rotation space, computing only the two 
leading LKCs, we can avoid these. The problem is not that the dimension of 
rotation space is higher, but that the Riemannian metric it induces is much 
more complicated than in the scale space case. Nevertheless, in Section 6 we 
shall see how once the geometry of the rotation space is well understood, 
the two leading Lipschitz-Killing curvatures of the scale space are effortlessly 
computed. 

Throughout, we have tried to carry out computations as carefully as pos- 
sible, so that not only can the interested reader follow the calculations, but 
can also see how they need to be adjusted for other cases. 

5. Proof of Theorem 2.1. We shall break the proof of Theorem 2.1 
into stages, the first of which relates to the crucial orthogonality properties 
mentioned in the Introduction. 

5.1. Orthogonality and the Riemannian metric on R N x {S}. Our aim 
in this subsection is to establish the orthogonality of the two subspaces 

T t R N , T 5 Sym+ C T t R N T s Sym+ ~ T [t>s) R N x Sym+ , 

under the inner product defined by the Riemannian metric induced by the 
process. Another factor that comes into play is a structured dependence of 
the restrictions of the Riemannian metric on R N x Synii to fibers of the 
form R N x {S}. This will turn out to be the key that allows us to carry out 
further detailed computations. 

To establish this independence, let (t, S) be a point in R N x Sym^- and 
let X, Y two derivations in T^R-^, T^Syin^ c Tr ti g\R N x Sym^ respectively. 
Using the connection to the covariance function given in (2.25), and writing 
Xt and YS to denote the elementwise derivatives of the identity functions 
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on R^ and Sym^ respectively, we have 

9( x,y) = ^ s) F (a) c(( t ,s),(t.s))| (is)Kfs) 

= 2\S\* (y\S\^ J ((XtySu)k(u'Su)k(u'Su)du 

+ 2\S\^ J ((XtySu) (u'(YS)u) [k(u'Su)] 2 du 

= 0, 

where the last equality follows since both integrands are odd functions. This 
establishes that the subspaces T t M. N and T^Syin^- are orthogonal at each 
point. 

We now compute the restriction of the Riemannian metric g to a fiber 
M. N x {S}, again based on (2.25). Writing, again, \?tf for the spatial deriva- 
tive of / with respect to the elements of t only, the Gram matrix of the 
Riemannian metric on R^ x {S}, relative to the standard Euclidean basis, 
is 



Var{V t /(i,S)} 

= \s\ 



4 \S\ 
4|S|*S 
4\S\^ S 

AS. 



V t k ((t - u)' S(t- u)) (V t k ((t - u)' S(t- u)))'du 
k (w'Sw) Sw (k (w'Sw) Sw\ dw 



\S[ 



k (w'Sw^J ww dw 
(k (v'v)) 



S 



2 -I , _i 
S 2 vv S zdv 



S 



k {y'v^ ) vv' 'dv 



Si 



5.2. The Riemannian metric on {t} x Sym^ . Relying on (2.25) to express 
the Riemannian metric and using a simple change of variables it can easily 
be seen that the restriction of g to {t} x Sym^ is independent of t. Denote 

this metric on Sym^ ~ {t} x Sym^ by g^ym^ _ Using the identity V5 \ S\ = 
\S\ S , valid on Sym^, we have 



V s (\S\*k((t-u)'S(t-i 

= \\S\*S- x k((t-u)'S(t-u)) 
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+ \S\H ((t - u)' S (t - u)) (t-u)(t-u)' . 
Therefore, for tangent vectors A,B£ Sym N ~ TsSym^, 

y s (\S\h((t-u)'S(t- u ))},A) 

V s (\S\H((t-u)' S(t-u))) ,B)du 



gf< (A, B) 



\S\ 



-k (u'Su) (S~\A) + k (u'Su) (uu',A) 
■k (u'Su) (5 -1 , B) + k (u'Su) (uu', B) 



du. 



After some slightly tedious calculus, and based on the fact that AC4 — 3«2,2 = 
(see (2.13)), all of the terms above can be expressed in terms of moments of 
the kernel k and its derivative k as follows. 
Sym^ 



9s 



(A,B) 



= Qj + § + «2,a) Tr(5- 1 ,4)Tr(5- 1 S) + 2 K2i2 Tr(S- 1 AS^S). 

When /i is the Gaussian kernel h(u) = e - ""'^/ 2 /ir N ^, and therefore 
k(x) = e~ x l 2 /-k n I a , simple calculations show that ^ + 2 + K 2,2 = 0, leaving 
only the term 2k2,2^(S~ 1 AS~ l B) above. 

5.3. Putting everything together. Let l-if 1 *^ , denote the Riemannian 
measures on A; x {S} ~ A induced by g and "Hp Hausdorff measure; i.e. 
Riemannian measure under the standard Euclidean metric on T. A standard 
calculation gives 



dU\ 



A,x{S} 



dH 



(t) = x l/2 \s\ 



TtdiA, 



where l^l^g^ is defined to be the determinant of PSP' for some I x N 
matrix P the rows of which form an orthonormal basis of TfdiA[ (under the 
Euclidean inner product). The determinant, of course, is independent of the 
choice of basis. 

Working in a chart and relying on the orthogonality established above, it 
is easy to see that 

C k+ i {At x B k ) = n k +i(Ai x B k ) 



(5.1) 



f [ drlf' X{S} (t) dhf" (S) 
Jd k B k JdtAi 



\S\\ dlAl dK B k k (S)dHf(t). 



d t Ai Jd k B k 
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The equality above relies on the expression for the Radon-Nikodym deriva- 
tive above and a Fubini theorem. 

We recall now our assumption that B k is orthogonally invariant, i.e. S G 
B k => QSQ' G B k for any Q G O(iV). Thus the mapping 5 i-> QSQ' is a 
diffeomorphism from to itself. Moreover, for the covariance function of 
the rotation field, we have 

C((t, S), (t, §)) = C((t, QSQ'), (t, QSQ')). 

Therefore, S h-> QSQ' is also a Riemannian isometry. Thus the measure 
T~L k k , induced by the Riemannian metric, is also invariant under S i— y QSQ'. 
This implies that we can drop the dependence of the integral in (5.1) on 
the tangent space Ttd\A\. That is, we can replace it, for all t, with a fixed 
j-dimensional subspace of our choice, or a random one. Therefore 

n k+l (A t x B k ) = x^-HifrAt) [ MS)dn^ (s) 

(5.2) Jd k B k 

= X^HiidiA^F^Bk). 

To complete the proof of the theorem, one must repeat essentially the 
same calculation for the codimension part of A\ x B k . 

6. Proof of Corollary 2.4. In this section we perform the computa- 
tions necessary to evaluate the expression in Corollary 2.4. These calcula- 
tions are much simpler than those for rotation space. The set B(<r^" 2 , a^ 2 ) 
can be stratified as 

<9iID>(ovrVr 2 ) = { vI :<J2 2 <v< (?i 2 } 

Clearly, each stratum is orthogonally invariant. Consequently, Theorem 2.1 
implies we need only compute Fi j jv(9iB(o"^" ,<rj~ )), Fo ) jv(QoB(o - 2 , )) 
andFi.jv-i^iB^ 2 ,^ 2 )). 
Clearly, 



F 0)N (doB(a- 2 , a~ 2 )) = a^ N + a~ N 



while 



F ltN (d 1 lD)(a^,a^)= / ^ v N l l f(y)dv 

2 2 

-2 



2 
2 
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where f(y) is the appropriate density with respect to Lebesgue measure. 

Above, we expressed the necessary integrals in terms of the parameteri- 
zation of <9iB(<7^~ 2 , cxT/ 2 ) by the curve v i— > v ■ I. Hence 

/(") = y^rhij) 

= v~ l ^N 2 + § + «2, 2 ) + 2NK 2 ,2 

Therefore, 

^-.(WVr 2 ) = c k , NJW ^ (*[ N - 1} - ^ 2 - (7V - 1) ) . 

7. Proof of Theorem 2.2. In this section we carry out the compu- 
tations necessary for evaluating the expressions in Theorem 2.2. Firstly, we 
describe a stratification of S(cr^~ 2 , cj -2 ), which is one of the regularity condi- 
tions required for applying the GKF. We shall actually only use the strata 
of codimension and 1 in our computations. However, if one wanted to com- 
pute higher order Lipschitz-Killing curvatures it would be necessary to use 
the full stratification. 

7.1. Stratification of the space. Our first step is to express the manifold 
M = T x §(<7^~ 2 , <T]~ 2 ) as a Whitney stratified manifold. Finding a Whitney 
stratification for ^(a^ 2 ,cr^ 2 ) itself yields a stratification for M as the prod- 
uct stratification, given that a Whitney stratification of T is assumed given 
as part of the data of the problem. Thus our first goal is to stratify S(ti, T2) 
for arbitrary < T\ < r 2 < 00. 

Semialgebraic orbits of the Lie group action GLn (R) x Sym^y — > Syrn^, 
(P, S) PSP', 

(7.1) {5 G Sym^vl rank(S) = Jfe, S>0},0<A:<iV 

constitute a Whitney stratification of {S G Sym^l S 1 > 0} (cf. [7], Page 21). 
Since 

,4 fe ( ri ) = {5 £ Sym^v | rank (S - tJ) = N - k, r x I < S} 
B k (r 2 ) = {S £ Sym N | rank (S - r 2 I) = N - k, S < t 2 I} 
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are obtained by applying affine transformations to the strata of (7.1), they 
form Whitney stratifications of 

(7.2) {S G Sym N | Tl I < S} and {S G Sym^ | S < t 2 I} . 

The final step in stratifying § (n, 75 ) is to verify that these stratifications 
intersect transversely. In order to do so, it is sufficient to show that for any 
j, k > such that j + k < N , & matrix S G A^(t\) n B k {r 2 ) and nonzero 
matrix C G Sym^y, there exist smooth paths A(t) : (—5,5) — > A 3 (t\), 
B(t) : (-5,5) — ► B k (T 2 ) with A (0) = B (0) = S such that for any t G 
M,5), 

(7.3) A (t) - 5 (t) = (A (t) - S)-(B (t) -S)= tC. 

Based on the fact that the group action mentioned above preserves rank, 
the problem of finding the paths A (t) , B (t) is equivalent to defining them 
by 

(7.4) A(t) = X(t)(S-nI)(X(t))' + n l, 

(7.5) B(t) = Y (t) (S - t 2 I) (Y (t))' + r 2 I, 

and finding appropriate invertible matrices X (t) , Y (t). 

Even under the restriction X (t) = Y (t) a solution exists. Substitution in 
(7.3) gives 

x (t) (s - n i) (x (t))' + Tl i-x (t) (s - t 2 i) (x (t)) 1 - t 2 i = tc, 

which is equivalent to 

X(t) (x(t))' = —!—C + I. 

t 2 - n 

For small t, define X (t) to be the unique positive definite square root of 
T jL Tl C + /. Substituting X (t) in (7.4) and (7.5) defines, for small 5, the 
required paths A(t), B (t), and proves transversality. 

The transversality of the intersections, together with the stratifications in 
(7.2), yield the Whitney stratification 

S (71,75)= |J (71,T 2 ), 

0<j,k 
j + k<N 

where §j> k (71,72) = A? (n) D B k (t 2 ). That is, we have found that the par- 
tition of S (71,72) according to the geometrical multiplicity of t\ and T2 is a 
Whitney stratification. 
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Note for later use that the submanifold of Sym^y of matrices of rank 
N - k, and thus A k (n) and £ fc (r 2 ), are of codimension k (k + 1) /2. Since, 
for j, k > with j + k < N, A^(t\) and B k {j2) intersect tranversely, the 
codimension of S j ' k (n, r 2 ) is [j (j + 1) + k (k + l)]/2. 

7.2. Volume density on Syni^. Having computed the Riemannian metric 
on Sym^, in order to apply Theorem 2.1 to the rotation space random field 

on T x §(0-2" ,<7^ ) we first explicitly compute ctH^ 2 ' Ul /dH^, i.e. the 
density with respect to Lebesgue measure on Sym^ when viewed as M K via 
the choice of a basis for the Hilbert-Schmidt inner product on Symjy. 

Before turning to the calculation we need some notation. Denote the stan- 
dard basis of Symjy by e = {eij)^l^ - =i . Normalizing the elements ejj with 

i 7^ j by a factor of 1/V2 yields an orthonormal basis relative to the Hilbert- 
Schmidt inner product which we denote by e = (e#)£^- =i . 

Working with the chart defined by mapping a matrix to its coordinates 
relative to e we see that the density is given by 

dn K ff s\_. A^l SymJ 



( 7 - 6 ) aoJe (S) = Wdet g^* (e n 

un-K V ^ J m=l,n=l 

where (e m )^ =1 = (e i ( m ) J -( m ))^ =1 is obtained by vectorizing (%)^^- =J 



As already stated in Subsection 5.3, for any Q 6 O (N), the mapping 5 1— > 

QSQ' is a Riemannian isometry of g^ mN , and thus preserves the measure 

' . Clearly, it also preserves the Euclidean Hausdorff measure "Hr. 
Therefore the density in (7.6) is invariant under S 1— >■ QSQ' as well. Since 
any real symmetric matrix is orthogonally diagonalizable, it follows that the 
density only depends on the eigenvalues of S. It is therefore sufficient to 
compute it for diagonal matrices. 

Let A = diag (Ai, . . . , \n)- For any two basis elements of e, 

Symt ,~ ~ , 
9 a {eij,e k i) 



+ k 2 ,2 + \p\ Tr (A- 1 ^) Tr (A" 1 ^) 
+2k 2i2 Tt (A-^j-A- 1 ^) 
^— + k 2 ,2 + -pj K 1 ^ 1 ^^^ + 2ft2,2A 4 ~ 1 AT 1 (5 ifc (5 77 . 



Suppose we order the basis in such a way that its first N elements are 
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ex i, . . . , ej\r at- Then, the corresponding K x K matrix has the form 



Sym+ 
9a 



K ' K (( i i 1* 



m=l,n=l 



where 

F= [{at^.o^.o 

V K—N 

By inspection, this is a rank-one perturbation of a diagonal matrix. Now, 
for a diagonal matrix D and vector u, det (D + un') = det (D) (l + v! D~ l u) , 
which, when applied in this setting yields 

-2 -2\ 



- (A) = Jdetj^^e;)}^ = , 
"'i/i' V ^ J m=l,n=l 

where q^jv is defined in (2.16). By the argument above, for general S € 

§KVr 2 ), 



(5) = c fe)JV |S| — 



7.3. T/ie push- forward to eigenvalues. Rewriting the integral defining 

Fk,n{^(^2 i a i ))' (2-26), in terms of the density dM.^" 1 ,<Tl /d%^ gives a 
relatively complicated integral. Fortunately, there is a simple way to convert 
it to an integral against a simpler density and, more importantly, over a 
much simpler domain. To do this, we need first to recall a classical result 
from the theory of random matrices. 

The starting point of the computation of the eigenvalues density in the 
Gaussian orthogonal ensemble (GOE) model (cf. [5]) is similar to our situ- 
ation. We begin with a measure fi on Sym^y with density (relative to 
that depends strictly on the eigenvalues. The goal is to compute the density 
of the push- forward of /U, denoted by u, under the mapping a : Sym^y — > K w 
which maps a matrix to the vector composed of its ordered eigenvalues. 

In the case of the GOE the density of fi is simply the product of the 
densities of the independent entries of the (real symmetric) random matrix. 
The entries of the matrix are zero-mean Gaussian variables with variance 2 
on the diagonal and variance 1 off the diagonal. That is, 

* (S) = 2~-f2-f(2,)-^e* P {-!y>- J £ S,} 

att K j=i l<i<j<N 
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n-v , , N(N+1) 

= 2—{2k) — exp{-±Tr(S 2 )}, 

where the factor of 2 N ~^ comes from the fact that was defined via 
idntification with the Hilbert-Schmidt basis and not the standard basis of 
symmetric matrices. 

The classical result is that 

| (A) _ (*)-* 2^ f[ f^S-S*. IA (A)| e Xp { - jf) A?}, 

where Z is Lebesgue measure on ~R. N and A (A) = Y\i<j (-\? ~ ^i) ^ s the Van- 
dermonde determinant. 

The specific Gaussian nature of fi has only a limited role to play in this 
result, and all that is really relevant for our needs is that the density is 
determined by a (S) at each point. In fact, for any measure Jl with such 
a density, we can write dfl/dM^ (S) = (^(a(S)) and rely on the GOE 
argument to conclude that for the corresponding push-forward measure v 
we have 



~" IT ^7§1a,<...<a„ I A (A)| Cf (A) , 



_ -'«/»> 

where fx and v are the measures corresponding to the GOE model. 

In particular, if we denote the push-forward of % K 2,1 by a with p, 
then 

N 



% (A) . (*r)*-»^n (^f^) ^ ^ IA (A)l 



C fc , Jv l Al <...<A JV |A(A)| J]A 



N N+l 
2 



3 



7.4. T/ie too leading LKC's. In view of Theorem 2.1, having computed 
the density -S (A), the computation of Cd (M) simply becomes a matter 
of rewriting the integral of Fk,n (pK& { a 2 2 ' a i 2 )) m terms of eigenvalues. 
Doing so leads to 



-2 -2^ 
2 >°1 y 



2 '"1 
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=c k , N / 2 \A(\)\n\]*d\. 

J<J 2 <Xl<---<^N<°' 1 j = l 

We now turn to F^ i 7v-i(<9^S((T^ 2 , 2 )), for which an identical argument 
yields 



N 



= C KN I I A (A) | 7^-i(diag(A)) J} \ ^ dX. 

Ja 2 <\i<---<\ N <a 1 - =1 

It remains to evaluate i ? ft-_i ] Ar(S((T^ 2 , <rj~ 2 )). Although this depends on 

rl K _ x , we shall not have to study this measure directly m order 

to compute the integral. Rather, an application of Federer's coarea formula 
([3, 6]) with a carefully chosen map will suffice. 

For any < r < a^ 2 define §mm( T > ^l" 2 ) *° ^ e the se ^ °f rea ^' symmetric 
matrices with unique minimal eigenvalue and with all eigenvalues in (r, crj~ 2 ) . 
Consider the smooth mapping \ m in '■ ^(O, c^ 2 ) — > M which maps a matrix 
to its minimal eigenvalue. To apply the coarea formula we need a Riemannian 

structure on the domain and range of X m i n - For the domain we take g^ mN 
as the Riemannian metric, and on M consider the Euclidean metric. 

Now, fix some < r < erf 2 . Let S G S^Jt, crj~ 2 ) and write S = QAQ' 
with Q £ O (N), A = diag(A) diagonal and A = (Ai, . . . , Aj\r) ordered so that 
Ai < ... < X N . 

Defining = Qe^ Q', e9 = ( e9- ) forms a basis of Sym^-. Therefore 

we can choose an orthonormal basis (bij)^l^ ' - =i of Ts§ min (t, erf 2 ), relative 
to g^ mN , such that 



e 3 Aie° 
e i,i A i e i,i 



where -D^.tv is defined in (2.17), and such that all bij with ^ (1, 1) are 
linear combinations of . . \ j ■ Since the basis is orthonormal, 

the Jacobian can be expressed by 



N,N 



(VA min (S) ,bij) 2 
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Note that if we can show that (v\ m in (S) is zero for any 7^ 
(1,1), then it will follow that 

JX m in {S) = |(VA m j n (S) , 611) I . 

In order to do so, and also compute this value, we examine the eigenvalues 
of S + tefp for small t. Write 

5 + tgg = Q(A + %)Q', 

and note that therefore we can work with the eigenvalues of Ut — A + teij 
instead. When i = j, Ut is diagonal with eigenvalues 

{Ai, . . . , Aj_i, Aj + t, Aj+i, . . . , Xn} ■ 

When i 7^ j, Ut also takes a very simple form and by a simple computation 
its eigenvalues are 

{Ai, . . . , Aj_i, Aj+i, . . . , Xj-i, Aj + i, . . . , An} 

U { \ {K + Xj ± ^{Xi-Xjf + 2t^j j . 

Therefore 

Xmin [S + teZ] - A mm (S) 
(VA mm (5) ,eg) = lim 5> ^ = 5(1,1) (*,i) , 

which implies 

<J "min [D ) — — . 

L>k,N 

Finally, Federer's coarea formula gives 

^(w 2 ) 



/ JX min (S)\S\*dH£ 
Js u . fr,o-r 2 ) 

Tim \ ' 1 / 



(S) 



= f dx [ \S\^ d7&"l^ 2 hs) . 

Since the integral on X~l in (x) is continuous in x, differentiating by r yields 
! i 2 |5|5 dH^ (r ' CTl " 2) (5) 

*' A mi«(°'2 2 ) 
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d_ 



[ j\ mm (s)\s\Un s i nm{T ' a ' 2) (s) 

=^ 2 J ^ in (r,^ 2 ) 



Ck. 



N 



N 



o-- 2 D k N Ja^<X 1 <---<\ N <a^ 



|A(A)|A 1 []A J 2 d\ 



3=1 



Ck,N N-2 

— °2 

Dk,N 



N-l 



A (A) | ] [ U; 2 lAj-ffj 



dA, 



' CT 2 — ^1— "' — -^iV-l^ "! j = l 

where the second equality is implied by the fact that any point of 

{ A : r < Ai < • • • < X N < af) \ a (S^ (r, a^ 2 )) 

is a boundary point or has zero density and the last equality follows by 
writing the integral as an iterated integral and differentiating. 

Obviously, we can also consider the maximal eigenvalue, define X max and 
^>max ( <T 2~ 2 ' r ) accordingly, and carry out a similar computation. Combined 
with the last result and the fact that dx-i^ (c^~ 2 , &\ 2 ) = A~^ n (c^~ 2 ) U 
Knax (°T 2 ) we then have 

Fk-i,n (9k-i§> (<?2 2 i a i 2 )) 

is^dnff'^is) 



9if-lS(o" 2 2 ,cr 1 2 ) 

Ck ' N f iA(A)inA;- 

k,N ^cr~ 2 <Ai<---<A A r_i<(T 1 7 2 - =1 

JV-1 N-l 
„N-2 IT I \ ^-21 , -JV-2 TT |„- 



AT-l 



3=1 



3=1 



dX. 



This completes the proof of Theorem 2.2. 



8. Scale space fields. In this section we shall compute all the Lipschitz- 
Killing curvatures for the scale space random field, under the assumptions 
of Section 2.2 on the kernel h and of Section 2.4 on the parameter space T. 
In this section we require only that h satisfy (2.7), without requiring full 
rotational symmetry. 

In the final analysis, the formulae that we shall obtain are quite simple. 
They involve no more than the Euclidean Lipschitz-Killing curvatures of T, 
the parameters A and k of (2.7) and (2.8), and combinatorial coefficients. 
On the other hand, there are many elements to each formula, so that any 
explicit computation of them will involve a computer. Thus, we shall proceed 
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by setting out two general formulae - (8.1) and (8.2) - and then in Section 
8.7 list explicit expressions for each of the summands in these formulae. 
Putting it all together to get actual numbers is a simple computing exercise. 

As far as the two leading Lipschitz-Killing curvatures are concerned, it 
requires no more than a few lines of algebra to see that the results of this 
section coincide with the earlier computation in Corollary 2.4 

For notational convenience in what follows, we make the simple scale 
transformation s = —In a, so that our random field is now denoted by 
f(s,t) and has covariance function 



The parameter space of /, which we shall denote by M, is now of the form 
T x [si, S2] C I^ 1 , for some —00 < si < s 2 < 00. 

8.1. Stratifying the parameter space. We start by stratifying the param- 
eter space M into manifolds of common dimension. In particular, we write 
the (N + l)-dimensional parameter space 



as the disjoint union of four types of pieces, arranged in three strata, ac- 
cording to their dimension: 



8.2. What we need to compute. Ultimately, we need to compute the 
Lipschitz-Killing curvatures Ci(M) of M under the Riemannian metric in- 
duced by the random field /. Rather than doing this directly, we shall com- 
pute certain other curvatures, denoted by jCf(M), specifically designed to 
be simpler to handle on sets of constant curvature. (We shall soon see that 
M has constant curvature, given by — k _1 .) 

There are simple relationships between the Ci(M) and £f(M), including 
the following one, which is (10.5.12) of [3]. 




M = T x [si,s 2 ] 



d N -iM 



d N +iM 
d N M 



A 



M° = ( Sl ,s 2 )xT°, 

( Sl ,s 2 )xdT U {s 2 }xT° U {si}xT° 

"side" U "top" U "bottom", 
{si} x 3T U {s 2 } x &T. 




N + l-i 
2 



(8.1) 
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To compute the Cf(M) themselves, we first write 

N+l 

(8.2) Cf(M)= £i(M;djM), 

j=N-l 

and then restrict ourselves to the case a = — Applying (10.7. 10) 3 of [3] 
along with the (yet to be proven) fact that M has constant curvature — 
gives that, for i < j, 

(8.3) C-^'iM-djM) 

Here, for each t £ djM, Zjj is a normally distributed random vector of 
dimension N + 1 — j in the space TtdjM^-, the orthogonal complement 
of TfdjM in R w+1 , S is the Riemannian shape operator of M, NfM is the 
normal cone to M at i, and Hj is the volume measure on djM corresponding 
to the Riemannian metric induced by the random field /. 
The next step is to compute each element in (8.3). 

8.3. The induced Riemannian metric. The next step is to identify the 
Riemannian metric g that the random field induces on M. In order to de- 
scribe this, however, we need first to choose a family of vector fields gener- 
ating the tangent bundle of M. 

We do this sequentially, starting with the set T, ignoring for the moment 
the scale component of the parameter space. For T, let 77 be the (Euclidean) 
outward unit normal vector field on dT, and extend this to a full (Euclidean) 
C 2 , orthonormal tangent vector bundle X\, . . . ,Xjy, on T, with Xjy = r\. 
Enlarge this to a vector bundle on all of M by adding the vector field v = 
d/ds, the field of tangent vectors in the scale direction. 

The Riemannian metric g at points (s, t) £ M can now be calculated on 
pairs of vectors from the above vector field using the variances (2.8) and 
(2.9) and the independence discussed there. These yield 

( 8 -4) 9(s,t)( x i,(s,t), x j,(s,t)) = Ae 2s (5jj, 

(8-5) 9{s,t){Xi^ s j),V( s j)) = 0, 

(8-6) g{s,t){ v {s,t),V{s,t)) = K , 



3 There is a small, but, for us, rather significant typo in (10.7.10) of [3], in that there is 
a minus sign missing before the re in re/ 2 /2. 
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where Sij is the Kronecker delta. We note, once again, the fact that (8.5) 
implies that the metric g is actually of product form is essential to all that 
follows. 

Note that the above three equations also imply that the stucture of the 
normal cones to M is well described by the initial, Euclidean, choice of 
vector fields. There is, of course, no normal cone at points t 6 dpf+iM, 
the interior of M. For the sets in d^M, one of the vector fields rj and v 
describes the normal geometry. In particular, along the side the normal is 
an s-dependent multiple of r\ and along the top and bottom the normals are 
constant multiples of v. 

Along ayv_i all normals are linear combinations of elements in the vector 
fields 7] and v. 

8.4. The Levi-Civita connection. The next step involves computing the 
second fundamental forms in (8.3) requires identifying the Levi-Civita con- 
nections. If we think of M as being embedded in a smooth open subset of 
]R Ar + 1 i with the same induced metric, then we shall write V for the Levi- 
Civita connection on M and V for the Levi-Civita connection on ambient 
set. 

For our purposes, we need only to know how Vt operates on vectors 
normal to M. Furthermore, since the Lie bracket 

(8.7) [Xi, v\ = 0, for all 1 < i < N, 

we need only compute V Xi Xj, V u Xi and V v v. 

To start, note that a straightforward application of Koszul's formula 4 and 
(8.4) and (8.5) imply, for i,j,k = l,...,N, that 

g(V x X 3 ,X k ) = \e 2s (V Xi X 1 ,X k }=g(V Xl X j ,X k ), 

where V is the standard Euclidean connection, and 

(8.8) g{V Xi X h v) = ~§-/ S &i3 = -Ae 2s <%. 

Now note that for any tangent vectors X, Y to T, Vx^ is a vector in l^" 1 " 1 , 
and so can be written as Ylf=i a j-^j + ^ u f° r appropriate coefficients. This 
fact, together with the last two equalities and (8.6), imply that for any two 
vector fields X, Y on T, 

(8.9) V X Y = V X Y - KT l Xe 2s {X, Y}u, 

4 This is (7.3.12) of [3], and in our case simplifies to 2g(V X Y, Z) = Xg(Y,Z) + 
Yg(X,Z)-Zg(X,Y). 
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giving us the first computation. 

As far as Vj,f is concerned, note first that another easy consequence of 
Koszul's formula is that 

(8.10) 0(V„i/,i/) = 1 -v{g{y % v)) = ^ 

giving us V^z/ = 0. 

All that remains is to compute V^Xj. An application of the Weingarten 
equation 5 and (8.9) yield 

g{V v X it X-) = | |;e 2s <% = \e 2s 5 i3 = g{Xi,Xf) 

and 

g{V v Xi,v) =0. 

Applying now the torsion freeness of connections 6 , along with (8.7) to the 
above, gives 

(8.11) V u X i= V X M = Xi, 

and we have the last of the three cases we were seeking. 

8.5. The second fundamental forms and the curvature matrix. With the 
relevant connections determined, we can now finally turn to computing sec- 
ond fundamental forms and curvature matrices along the three stratifications 
of M. 

1. 8n+\M, the interior (si,s<2) x T° : 

Since the normal space in the interior is empty, there are no normals, no 
second fundamental form and no curvature matrix here. 

2. d N M: The side (si,s 2 ) x dT: 

A convenient choice for an orthonormal (in the metric g) basis for the 
tangent space at any point on the side is given by 

f x l X N ^ v x 
lAVV''"' AVV' kV2 J ' 

with outward unit normal vector rj/\ l / 2 e s . 

5 This is (7.5.12) of [3] and states that the scalar second fundamental form is given by 
Su{X,Y) = 5 (Vxi» = -g(Y,Vxv). 

6 This is (7.3.10) of [3], and states that V X Y - V Y X - [X, Y] = 0. 
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The scalar second fundamental forms of interest are therefore 



all of which can be computed directly using Weingarten's equation and (8.9)- 

(8.11) . We summarise them in a curvature matrix, writing 

C = \- 1 / 2 e- s C t 


where Ct is the (N — 1) x (N — 1) Euclidean curvature matrix of dT in the 
basis {Xi, . . . , Xtv-i, rf\. Recall for later use that, from Section 7.2 of [3], in 
the current scenario, the trace in (8.3) can be replaced by 

(8.12) (j - iy.detvj-iC, 

where detrj is our usual sum of determinants of j x j principal minors. 

3. d N M, the bottom {si} x T° and the top {s 2 } x T° : 

Starting with the bottom {s\} x T°, the outward unit normal vector is 
and a convenient orthonormal basis is given by 

X\ Xn 



X l/2 e s ' ' ' ' ' A l/2g S 

The curvature matrix is therefore N x N with entries 

by Weingarten's equation and (8.8). 

Thus, for the bottom, the outward curvature matrix is k~ 1 I 2 InxN while 
along the top the same arguments give it as — k~ 1 I 2 In x n . 

4. d N ^iM: The edges {si} x dT and {s 2 } x dT: 

For the edges, we need to consider the scalar second fundamental form 
itself, and not just the curvature matrix. As above, an orthonormal basis for 
the tangent space is 

X\ Xn_i i 



AVV''"' AVV' «V2_ 
but now an orthonormal basis for the normal space is {X~ 1 ^ 2 e~ s n, k~ 1 I 2 v}. 
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Applying the Weingarten equation, we need to compute 

9 (v A -i/ 2e - SXi (X-^e-'Xj) , aA" 1 / V s r? + 6k" 1 / 2 / 
for arbitrary a, b. Applying (8.9) gives that this is 

where, with a minor abuse of notation, we now use Ct to denote the (N — 
1) x (iV — 1) Euclidean curvature matrix of dM at t. 

8.6. The curvature tensor. We now have all the pieces we need to com- 
pute the Lipschitz-Killing curvatures C~ K , and could actually proceed to 
the final computation. However, in justifying the formula (8.3) we used the 
fact that M has constant negative curvature — k . Now we shall take a 
moment to prove this. 

Let 

R(X, Y) = V x Vy - Vy V x - V[X,Y] 

be the curvature operator. Using our previous calculations of the connection 
we have 

R(Xi, v)X k = V Xt V u X k - V u V x X k 

= V x X k - V v {V x X k - SikK^Xe^u) 

d 

= V x X k - 5 lk K~ l Xe 2s v - V x X k + 5 ik K- 1 — Xe 2s u 
= d ik K Xe v, 

R(Xi, v)v = V Xi V v u - V v V Xi u 
= -VuX l 
= —Xi, 

R(Xi,Xj)u = V Xi V Xj u - V Xj V Xi u - V [x . iXj] u 
= V x X J -V x X t -[X l ,X j ] 
= V x Xj - V Xj Xi - [X it Xj] - K^Xe^iSij - b 3i )v 
= 0, 
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R(Xi,Xj)X k = Vx^XjXk - VxjVxi^fc - V[Xi,Xj]Xk 

= V Xi (V Xj X k - 5 3k K- l \e 2s u) - V X] (V Xi X k - S^Xe^u) 

- V [Xi ,x s \Xk + ^Xe^dX^X^X^u 

= V Xi V Xj X k - 5 jk K~ l Xe 2s Xi - K^Xe^iVxjX^X^u 

- V Xj V Xi X k + 5 lkK ~ 1 Xe 2s X j - K^Xe^iVx^X^u 

- V { x irXj \X k + K-WdX^X&Xju 
= VxiV^Xfe — VxjVxiX k - V[x u Xj\X k 

+ ({[XuX^Xt) - (V Xj X k ,Xi) + {V Xi X k ,Xj)) 

- 5 jk n- l Xe 2s Xi + SikK^Xe^Xj 

= -5 jk K~ 1 Xe 2s X i + 5 ik K~ 1 Xe 2s X j . 

where the last equality follows from the flatness of Euclidean space and 
torsion freeness. 

Since the curvature tensor is given by 

R(X,Y,Z,W) = g(R(X,Y)Z,W), 

it is now easy to use the above calculations to check cases and see that 

R{X, Y, Z, W) = -k~ x I 2 {X, Y, Z, W)/2, 

where / is the identity form given by I(X, Y) = g(X, Y), and, with the usual 

tensor product I 2 (X, Y, Z, W) = I(X, Z)I(Y, W) - I{X, W)I(Y, Z). 

From this, of course, follows our claim that M is a space of constant 
curvature — k~ 1 . 

8.7. The Lipschitz- Killing curvatures Cr K 1 . With all the preparation 
done, we can now begin the computation of the Lipschitz-Killing curvatures 
Cj K (M; d k M), according to their definition in (8.3). As in the discussion of 
second fundamental forms, we divide the computation into separate sections, 
each corresponding to a different stratum in the stratification of M. 

Throughout the following computation we take Z\ and Z2 to be two 
independent N(0, 1) random variables. We shall associate the Zj with the 
normal vector fields to obtain normal vectors of the form Z\v + Z217. 

We also adopt the notation 



5.13) n k = E{z}l Zj >o} 



(\ if k = 

2n "^l if k = 2n + 1 is odd 

Z7T 



(2n-l)(2n-3) 



if k = 2n is even. 
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1. dN+\M, the interior (s\,S2) x T° : 

Since the second fundamental form is zero in 8m+\M, the only non-zero 
Lipschitz-Killing curvature occurs when N + 1 = j = i in (8.3). In this case 
we have 

£j£i' (M; a N+1 M) = «,„((«, s,)xr] 

(8.14) = " U2Xm J" l T /' dtd ' 

Ns 2 _ N Sl 

= K i/2 X N/2 e ^ £ E N (T;T°), 

where C^(T;T°) is computed in the standard Euclidean sense. 

2. d N M: The side (si,s 2 ) x dT: 

For this case, all the Lipschitz-Killing curvatures need to be computed, 
and we replace the trace in (8.3) by the determinants of the curvature matrix 
as in (8.12). Then 

(8.15) £j K ~\M;( Sl ,s 2 )xdT) 

= (2vr)-(^)/ 2 E{Z 2 ^l {Z2>0} } 

f , (e~ s C t 0\ „, , \ 

x / detrjv-,- n ) dn N {s,t) 

J( S1 ,s 2 )xdT V u u / 

x f 2 [ (X-^e-^-^X^e'^detTN^iCt) dtds 

Js L JdT 

„(j-l)s2 — l)si 

= K V2 A (i-i)/2 « U > -e* > £ e ^ 

j - 1 J 

in a parallel notation to (8.14). 

3. d N M, the bottom {si} x T° and the top {s 2 } x T°: 
Beginning with the bottom, {si} x T°, for < j < N we have 

(8.16) £j" _1 (M;{ ai } xP) 

= (2 7r )-^')/ 2 E{^l {Xl>0} } 

x / detr7v-j(K 1/2 /7Vx7v)^Jv(t) 

J{si}xT° 

= (2^- 1 )-( Jv ^)/ 2 m - J A^e^^f (T; T°). 
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A similar result holds for the top, {^2} x T°, viz. 
(8.17) Cj^ 1 (M;{s 2 } xf) 

= (-l)^(2vr K - 1 )-^)/^^_ J ^A 7V / 2 e 7Vs2 £^(T;r ). 

4. d N -!M: The edges {si} x dT and {s 2 } x dT: 
We start with the top edge, {s 2 } x dT . 

£j K ' 1 (M;{s 2 }xdT) 

= (2vr)-( 7V - 1 -^/ 2 I Ejdetr^.L^Z^- 1 / 2 / - Z^T^e'C^ 

J{s 2 }xdT 

xl {z 1 >o} 1 {z 2 >o}}d'HN-i(s,t). 
Use now the easily checked expansion that, for < k < n, 

k , N 

detr fc (al nxn + A nxn ) = ^ a k ~ m ( _ J detr m (A), 

m=0 ^ m ^ 

to expand the detr term in the expectation above and see that 
(8.18yCj K ~ 1 (M;{ S2 }xdT) 

N ^Z j / N 1 m 

= ^y{N-X-j)/2 K -{N-\-j-m)/2l 



x 



ft 

\ 7 

m=0 v J 

E«- 1 ^-l { ^ 1>0} }(A 1 / 2 e-) JV - 1 E{^l { ^ 2> o } } 
x / detr m (A- 1 / 2 e- S2 C t ) dH N -\(t) 

JdT 



N-l-j 



(2vr K - 1 )-( 7V - 1 -^ 2 £ (AV2 e s 2) JV-i-m fN-l-m\ m/2 

m=0 \ 3 J 

x/ijv-i-j-^Z^l/^^o}} / detr m (Ci) dH N -i{t) 

JdT 

N-l-j /AJ . 1 s 
K (iV-l-i)/2 (A l/2 eS2) AT-l £ P- 1_m j(A- 1 /2 e -'»« 1 /2) 

™— n \ 3 J 



m=0 

x^ A r_ 1 „ J _ m £^_ 1 _ m (T;OT). 

A similar argument also works for the bottom edge {si} x dT, the only 
change being that the condition l^ Zl>0 y becomes 1|^ 1<0 }, giving the final 
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form 



(8.19^7" \M;{ Sl }xdT) 



= K 




N - 1 - m 



J 




in 



m=0 



„ / i \N—\—j—m,. 
X(-l) /Xjv-l-j-j 




(T-dT). 



Collecting (8.14)-(8.19) now gives us all the £~ K (M; djM), from which, 
via (8.1) and (8.2), we can compute the Ci(M), as promised at the beginning 
of the section. 
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